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Ultrafast magnetization dynamics in metallic heterostructures consists of a combination of local 
demagnetization in the ferromagnetic constituent and spin-dependent transport contributions within 
and in between the constituents. Separation of these local and non-local contributions is essential 
to obtain microscopic understanding and for potential applications of the underlying microscopic 
processes. By comparing the ultrafast changes of the polarization rotation and cllipticity in the 
magneto-optical Kerr effect (MOKE) we observe a time-dependent magnetization profile M(z,t) 
in Co/Cu(001) films by exploiting the effective depth sensitivity of the method. By analyzing 
the spatio-temporal correlation of these profiles we find that on time scales before hot electron 
thermalization (< 100 fs) the transient magnetization of Co films is governed by spin-dependent 
transport effects, while after hot electron thermalization (> 200 fs) local spin-flip processes dominate. 

PACS numbers: 75.78.Jp, 72.25.Ba, 78.20.Ls 


I. INTRODUCTION 

Electronic excitations in ferromagnets are essentially 
spin polarized. In the transition metal ferromag¬ 
nets Fe, Co, and Ni excitations of 0.1-1 eV relax on 
femto- to picosecond time scales HI due to scattering 
with secondary excitations mediated by electron-phonon, 
electron-magnon, and electron-electron / exchange inter¬ 
action. The underlying microscopic processes are essen¬ 
tial in ultrafast magnetization dynamics 3,0. 

Recently this field has been propelled by ultrafast spin 
polarized ml 117)1 and unpolarized currents 0 generated 
by a gradient in excitation density 0• These effects offer 
to exploit non-local magnetization dynamics. Already 
now femtosecond (fs) laser-excited spin currents in layer 
stacks were reported to induce spin transfer torque [13|. 
Furthermore, laser-induced spin currents in a magnetic 
tunnel junction were controlled by a bias voltage 0 - 
Moreover, spin currents are reported to drive an ultrafast 
change between ferri- and ferromagnetic order 0. 

Such demonstrations suggest an extension of spin- 
tronics, including spin filter and magnetoresistance ef¬ 
fects, into the non-equilibrium regime. In turn, fs time- 
resolved experiments allow conclusions on spin transport 
to shed light on the underlying microscopic 
mechanisms. However, up to date experimental means 
to distinguish ultrafast spin currents, which propagate 
into conducting contacts/substrates from ferromagnetic 
layers, and competing depolarizing spin-flip scattering, 
which reduces the magnetization M in these spin cur¬ 
rent sources, are missing. 
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Here we show how to fill this gap by establishing that 
laser-excited spin currents result in a specific transient 
magnetization profile M(z,t) in the direction z along the 
film normal, before spin flips impact the dynamics. We 
employ ferromagnetic Co/Cu(001) films as a model sys¬ 
tem for a spin current source contacted by a spin sink. 
Key is the combined analysis of the film thickness depen¬ 
dence with the time-dependent difference in polarization 
rotation 6 and cllipticity e of the magneto-optical Kerr 
effect (MOKE), which provides a depth sensitivity and 
access to We separate spin currents generated 

by non-thermalized carriers up to ^100 fs and local de¬ 
magnetization due to spin-flip scattering of thermalized 
carriers dominant after ~200 fs. 

A contested issue in ultrafast MOKE is whether the 
observations reflect M(t ), or transient changes of the op¬ 
tical constants. Koopmans et al. 0], finding that 9 
and e of the complex Kerr angle $ = 9 + i ■ e showed 
a different behavior for t < 100 fs, assigned the differ¬ 
ence to state-filling effects. Later, Guidoni et al. [l6j 
reported that the magneto-optical response is dominated 
by the magnetization dynamics after electron thermal¬ 
ization. However, the origin of the transient difference 
between 9 and e before thermalization has so far been 
unclear 0] . We derive that the transient difference of 9 
and e on femtosecond time scales is a result of a spatial 
profile in the magnetization. 

Also, the relative importance of (i) spin-flip scatter¬ 
ing of thermalized electrons, described by the micro¬ 
scopic three temperature model (M3TM) 0 , and (ii) 
superdiffusive spin transport 3 0, 0 of hot, non¬ 
equilibrium electrons, in ultrafast magnetization dynam¬ 
ics is under discussion. Experiments on layered struc¬ 
tures have demonstrated the importance of spin trans¬ 
port 331 , 0 , while demagnetization of metallic films 
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on insulators was ascribed to spin flips [2Cj , leading to an 
emerging consensus that both processes play a role [2lj ■ 
Here, we investigate both spin-flip processes confined to 
the ferromagnetic film, as well as propagating spin cur¬ 
rents that are accepted by the metallic substrate, in a 
model system. 


II. ANALYSIS OF TRANSIENT 
MAGNETIZATION GRADIENTS 

Femtosecond laser excitation of a heterostructure like 
Co/Cu(001) epitaxial films, which present a large gradi¬ 
ent in spin polarization from the film to the substrate 
already in equilibrium, leads to spatial redistribution of 
charge carriers dHHH and spin polarization by transport 
effects in particular across the interface [y, 0,|2l|]. Locally 
within the ferromagnetic film the magnetization changes 
also by spin-flip scattering. In order to quantify the spin- 
dependent transport processes compared to local demag¬ 
netization we first estimate these two contributions by 
analytical calculations in one dimension along the inter¬ 
face normal direction. Lateral magnetization gradients 
within the film plane are also present due to the lateral 
gradient in excitation density following the laser focus 
with a full width at half maximum, in the present case, 
of about 16 pm. However, since these lateral gradients 
are three orders of magnitude smaller than the one along 
the normal direction, we neglect them here and consider 
only the normal component, which is essential for the 
sub-picosecond dynamics. 


A. Estimation of spin diffusion and local 
demagnetization contributions 


In Fig. [lja) we contrast the relative local magnetiza¬ 
tion m(z,t) = M(z,t)/M 0 for (i) spin-flip scattering of 
thermalized electrons and (ii) spin transport, calculated 
with an extended M3TM and a spin diffusion equation, 
respectively. Mq is the homogeneous equilibrium magne¬ 
tization. 

To estimate the spin-flip contribution we start with cal¬ 
culating the time- and spatially dependent electronic and 
lattice temperatures, T e (z,t ) and T\(z,t), respectively, 
for different Co film thicknesses d on Cu(001) employ¬ 
ing a two temperature model (2TM) [23] • We include 
electronic heat diffusion, which is driven by the spatial 
gradient in the electronic temperature [23} , and account 
for the Co film with thickness d. We further include the 
Cu substrate in the calculation. We take the optical pen¬ 
etration depth (5 s kin of the pump pulse at 800 nm central 
wave length into account and include thereby the position 
dependent optical absorption in the calculation. The ini¬ 
tial excitation profile in cobalt and copper was calculated 

by 


S(z,t) = 


Z7 1 i . . p 2 /^skin ,i 

■ L abs ,2 ° 

<5skin,i • (1 - 


G(t) (1) 




FIG. 1. (Top) Kerr rotation and ellipticity of Co/Cu(001) 
are measured in a pump-probe experiment, (a) Normalized 
magnetization profile m(z) calculated with the M3TM 0 
(black) for to = 170 fs and spin diffusion (SD, red) for to = 
20 fs. (b) MOKE sensitivity £(z) for the real part/rotation 
(g(z) (black) and imaginary part/ellipticity ( e (z) (red). The 
respective curves are shown for Co thicknesses d of 6 (dotted), 
10 (dashed), 20 nm (solid lines). 


G{t ) is a normalized Gaussian function with a FWHM 
of 35 fs, di is the film thickness, where i stands for 
cobalt and copper. For copper we used dc u = oo. 
The penetration depth of the laser <5 s kin,Co=13 nm and 
<5skin,Cu=12 nm was calculated according to the optical 
constants nc 0 = 2.56 + i ■ 4.92 and ncu = 0.26 + i ■ 5.26 
[H|. These optical constants were also used to calculate 
the absorbed fluence .F abSi j employing [25|, where we set 
the incident fluence to F ln = 3 mJ/cm 2 . 

Subsequently, we used the M3TM, which considers 
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spin-flip scattering only for a thermalized electronic sys¬ 
tem, to obtain g?tom 3 Tm/ dt, and tom 3 Tm(z, i) E2- To 

take the thermalization time of the electronic system in 
the simulation for the spin-flip scattering into account, 
we multiply the expression drriM 3 TM/dt from the M3TM 
na with a time dependent factor ntherm(i) [22]. It 
was calculated with n^ot/dt = G(t ) — nhot/Ttherm and 
dntherm/dt = n ho t/Ttherm, where n hot {t) and n t herm(t) are 
factors between 0 < n < 1, G(t) is a normalized Gaussian 
function with a FWHM of 35 fs and Ttherm = 150 fs is 
the thermalization time of the electronic system, which 
is a typical time for a 3d metal [2f|. The used param¬ 
eters are the electron-phonon coupling constants geo = 
9.3 • 10 7 W/m 3 • K, gcu = 10 7 W/m 3 • K [13, the elec¬ 
tronic heat capacity coefficients yco = 4.4 mJ/mol • K 2 , 
7Cu = 0.69 mJ/mol-K 2 [13, the specific heats for T —> oo 
Cco = 24.81 J/mol-K, Cqu = 24.43 J/mol-K, the Debye 
temperatures ©c 0 = 386 K, 0c u = 310 K f29|, the mass 
densities pco = 8.86 g/cm 3 , peu = 8.96 g/cm 3 jjjlj], the 
molar masses Mq 0 = 58.93 g/mol, Mq u = 63.55 g/mol 
[H|, the thermal conductivities k c 0 = 100 W/m • K, 
KCu = 400 W/m • K [2£j, the Curie temperature of Co 
Tq = 1388 K [Tjj. and the M3TM scaling factor for the 
demagnetization rate of Co R = 25.3 ps” 1 (l7j . 

To simulate spin diffusion (SD) we employ a diffusion 
equation for excited electrons based on Fick’s second law 
|3J] for majority and minority electrons. Eq. [2] con¬ 
tains the respective source, diffusion and decay terms. 
We discard charge current contributions because con¬ 
trary to dielectrics a potential charge current which is 
linked to the spin current is screened in metals already 
on time scales slower than the inverse plasma frequency. 
In Co these time scales were concluded to be sub-fs [30j 
and are thus well below our experimental time resolu¬ 
tion. We use the diffusion coefficient D a = v a \ a /?>, a 
represents majority and minority electrons, for self dif¬ 
fusion to simulate the transport. The respective mean 
free path along z is A CT = v a r a . Velocities v a and life¬ 
times Tfj were averaged up to 0.5 eV above the Fermi 
level Ep which is reasonable to account for spin diffusion 
considering that the experimentally determined scatter¬ 
ing times reported in Refs. HI refer to primary excited 
electrons. An electron at an energy of E — Ep = 1 eV 
has, e. g., a lifetime below 10 fs. Within our time resolu¬ 
tion, such a primary excited electron will have scattered 
more than once and lowered its energy accordingly. We 
used energy averaged lifetimes r-j- = 22 fs and tj. = 20 fs 
which were determined from Ref. !■ The velocities 
were calculated by density functional calculations (see 
Appendix) which resulted in energy averaged velocities 
Uf = 0.6 nm/fs and uj, = 0.2 nm/fs. Taking these values 
for A and v we obtained the considered diffusion constants 
D-f- = 2.4 m 2 /fs,D| = 0.27 m 2 /fs. We finally calculate 
the spin diffusion by 


dn a 

~dt 


dn 2 

= Da ' + Ca ' S a (z, t) 


rig 
Tel ’ 


( 2 ) 


where n a is the density of excited electrons, c a Sg(z,t) 
describes fs laser excited primary and subsequently gen¬ 
erated secondary electrons with c a being a constant fac¬ 
tor. To account for the spatial distribution, we convo¬ 
luted two exponential decay terms according to the po¬ 
sition dependent absorption due to the optical penetra¬ 
tion length S s kin and the spin dependent mean free path 
Abaii.c- of primary excited electrons [HJ. The latter is de¬ 
termined by ballistic velocities and lifetimes which were 
taken at energies of the maximum excitation probabil¬ 
ity in the joint density of states, see Appendix. The 
values are Ubaii.t = 1 nm/fs, Ubaii.i = 0.3 nm/fs and 
Tbaii.t = 18 fs, Tbai4 = 6 fs. As an effective decay time of 
the spin polarized current we took a value of r e \ = 100 fs 
[ill ]. Then, Abaii.tr was calculated with the respective ve¬ 
locities and lifetimes. We included propagation in Cu 
with A = 70 nm [36( by setting n a = 0 since bulk Cu acts 
as an efficient electron and spin drain. The last term 
—ITei describes electron thermalization and spin cur¬ 
rent decay by a time constant r e \. The change dmsp>/dt 
is defined through the balance of the magnetic moments 
of in- and outgoing electrons. 

Comparing the magnetization profiles resulting from 
local spin-flip scattering and spin transport, as depicted 
in Fig. □la), m(z,to) are strikingly and systematically 
different for the two scenarios because the respective gra¬ 
dients dm(z , to)/dz have an opposite sign. While consid¬ 
eration of of spin-flip scattering of thermalized electrons 
results in a minimum m at the surface and a maximum at 
the Co-Cu interface, the spin diffusion description leads 
to a depletion of m at the interface in combination with a 
weak variation near the surface, in good agreement with 
a more sophisticated description by superdiffusive spin 
transport fl8j . 


B. Measurement of the spatio-temporal 
magnetization dynamics 

To analyze m(z, t ) experimentally, we measured the 
complex MOKE = 9+i-e in a pump-probe experiment, 
see Fig.QJtop). The magneto-optical (MO) Kerr rotation 
6 and the ellipticity e were determined by a polarization 
analysis using a balance detection scheme. We employed 
a cavity dumped Ti:sapphire oscillator which generates 
p-polarized 35 fs laser pulses at hi/ = 1.55 eV and 40 nJ 
pulse energy, that are split into pump and probe pulses 
at a 4:1 ratio, at 2.53 MHz repetition rate. The incident 
pump fluence was 6 mJ/cm 2 . Further details of the ex¬ 
perimental setup are given in Ref. [§3, El ■ Epitaxial Co 
films of 4 nm < d < 20 nm, which we investigated in 
situ , were grown in ultrahigh vacuum on Cu(001) follow¬ 
ing Ref. [391 ]. The film’s easy axis of the magnetization 
lies in the film plane and we measure MOKE in the lon¬ 
gitudinal geometry. 

We find different transient responses of 9{t) and e(t), 
shown in Fig. [2] Our following conclusions are based on 
the fact that Kerr rotation 9 and ellipticity e exhibit dif- 
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FIG. 2. Time-dependent relative MO responses of Kerr ro¬ 
tation A6(t)/6o (dots) and ellipticity Ae{t)/eo (diamonds) of 
Co/Cu(001) for d = 6, 10, 20 nm. Solid lines represent analyt¬ 
ical fit functions, while dotted lines guide the eye. For deter¬ 
mination of t = 0 we measured the auto correlation function 
with second harmonic generation from the surface (dashed 
line, bottom panel). Arrows indicate the delay of minimum 
MO response; t* indicates when A9(t)/9o = Ae(t)/eo- 


ferent depth sensitivities which we use to identify specific 
magnetization profiles m(z, f) caused by spin-flip scatter¬ 
ing of thermalized electrons or spin polarized currents, 
see Fig. [2(a). In Fig. [2b) we plot ( = (g + i ■ which 
we term MOKE sensitivity for the complex Kerr angle 
<J) = 6 + * • e, as a function of the position z within a 6, 
10, and 20 nm thick film. The sensitivity ( indicates how 
changes A m(z, t ) = m(z, t) — l contribute at a depth z in 
first order to the total change. The measured transient 
changes in $ are calculated by integration over the Co 
film thickness 


A9(t)/9 0 + i ■ Ae(f)/eo = f ( (z) ■ Am(z,t)dz. (3) 

Jo 


Following Traeger et al. [13, ((z) was calculated by 
setting m = 0 with the exception of a part z, z+dz where 
to = 1, and determining the complex Kerr angle of this 
part &dz = 9dz + itdz for this system. To take pump- 
induced changes into account we normalized $>dz to the 
equilibrium values and obtain 


C(z) • dz 


9dz 

ftT 


£dz 

£o 


(4) 


We used the matrix formalism by Zak et al. [4l[ with 
dz = 0.2 nm and refractive indices n c u and nco- The 
static MO constant qc 0 which enters the determination of 
C was determined experimentally to qc a = 0.017 —TO.020 
by fitting the thickness dependent MO contrast and the 
ratio e/9 for d = 2 — 20 nm (not shown). 

Note that </g is larger near the surface than near the 
Co-Cu interface and C, t exhibits a weaker z-dependence, 
see Fig. [2b). Applying eq. © reveals that (i) 9 probes 
effectively the near surface part. For thicker layers like 
20 nm the sign of the rotation sensitivity eliminates some 
signal contribution from the bulk part of the film which 
also leads to an effective probing of the surface near re¬ 
gion. (ii) The second magneto-optical observable e rather 
averages over the film and includes a sensitivity at the 
Co/Cu interface. This depth dependence of ( provides 
us with a probe for m(z, t). 

Now we consider our experimental results in detail. 
Fig. [2 depicts the time dependent relative MO re¬ 
sponses A9(t)/9o and Ae(f)/eo for Co/Cu(001) for d = 
6,10, 20 nm. The curves show a reduction of the signal 
which begins within the laser pulse and a subsequent re¬ 
covery starting between 100 and 300 fs, depending on 
the film thickness. Most importantly, we find for all d 
differences for 9{t) and e{t). For d > 10 nm we observe 
a faster reduction of e than for 9 and a crossing at delay 
t* of both curves within 200 fs. After this crossing the 
changes are inverted relative to each other and A9(t)/9o 
remains stronger than Ae(t)/eo up to 25 ps. This re¬ 
maining difference will be addressed further below. For 
d < 10 nm the initially stronger reduction of e compared 
to 9 is not observed, likely due to a limited pulse du¬ 
ration and a remainder of the coherent artifact (42j, in 
combination with a weaker dm/dz. Changes for such d 
are larger in 9 than in e. 

From the observed 9{t) and e(t) we conclude on m(z, t) 
and relate this to a spin polarized current or spin-flip 
scattering, see Fig. ©a). Here, the different depth sensi¬ 
tivity of 9 and e, see Fig. ©b), is essential. Remember, 
e exhibits a stronger sensitivity at the Co/Cu interface 
than 9. The stronger reduction of the ellipticity com¬ 
pared to the rotation starting within the pulse duration, 
see Fig. [2[ thus implies that the film is demagnetized 
more near the inner interface than near the surface. The 
respective magnetization profile agrees qualitatively with 
the one obtained for SD, see Fig. ©a). After some tens 
of femtoseconds, depending on d 1 m(z,t) is concluded 
to be redistributed towards a profile expected from the 
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FIG. 3. (a) Time delay r m i n of minimal response in A9(t)/9o 
(black dots) and Ae{t)/to (red diamonds) vs. d. Solid lines 
are linear fits, dashed lines depict simulations. 


M3TM p3, as can be seen in the larger variation of 9 
compared to e in Fig. [2] at t > t*, taking into account 
that 9 has a higher sensitivity at the surface than e, 
see Fig. mb). These transient changes of the magnetiza¬ 
tion profile suggest that when the electronic system has 
not yet thermalized, a spin polarized current dominates 
m(z,t). With electron thermalization, spin-flip events 
take over. This is consistent with both the superdiffusive 
transport model [18] and the M3TM [l7j|. The superdif¬ 
fusive transport model predicts spin transport for a non- 
thermalized electron system, which recedes with thermal¬ 
ization [43J, while the M3TM considers only thermalized 
electrons, which contribute to spin-flip scattering E3- 

Our conclusions are corroborated by finding that the 
time delay r m j n of the minimum MO response exhibits a 
systematic thickness dependence, as depicted in Fig. [2] 
To determine this delay of the minimum response we 
performed a fitting analysis and determined r m j n using 

/(f) = A m ■ (1 — exp(—t/rM)) + Ari ■ (1 — exp(—£/ tri)) + 

Ar 2 ■ (1 — exp(—t/r^)); A and r are the amplitudes 
and time constants, respectively, M,R1,R2 refer to de¬ 
magnetization and recovery [351 ]. For e and 9 the mini¬ 
mum positions shift to later times with increasing d, see 
Fig. [3] Particularly interesting is the different thickness 
dependence of e and 9. By fitting the minimum’s posi¬ 
tion as a function of d with a line, we find with slopes 
of bg = 13 fs/nm and b t = 4 fs/nm a pronounced, three 
time difference. 

Fig. H] shows time-dependent MOKE data for longer 
time delays up to 25 ps for selected him thicknesses. The 
difference between the MOKE ellipticity and rotation re¬ 
mains clearly longer than 5 ps and the time delay of 
merging of the two observables shifts to later delays with 
increasing him thickness. Considering typical electronic 
relaxation or electron-lattice equilibration times of few ps 



FIG. 4. Time-dependent relative magneto-optical responses 
of Kerr rotation A9(t)/9o (dots) and ellipticity Ae(t)/eo (di¬ 
amonds) of Co/Cu(001) for d = 6, 10, 20 nm for pump-probe 
delays of up to 25 ps. Solid lines represent analytical fits 
described in the text. 

suggests that this difference cannot be fully explained by 
the competition of spin diffusion and local demagnetiza¬ 
tion. We will come back to this aspect in the discussion 
section. 


C. Model description of the spatio-temporal 
magnetization dynamics 

We simulated the thickness dependent m(z , t) with the 
extended M3TM ’17] and the spin polarized diffusion 
equation, as introduced above, by adding both contri¬ 
butions 

dm(z,t) dm M 3TM{z,t ) dm S D(z,t) . . 
dt dt + dt ' U 

We calculated the corresponding MO response with 
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eq. fusing the parameters given above in Sec. II.A, and 
determined the delays of minimum MO response as a 
function of film thickness d, which are included in Fig. [3] 
as dashed lines. The overall shift of r m j n to later times 
with increasing d is explained bv the lower heat conduc¬ 
tivity of Co compared to Cu [23. Essential for the dis¬ 
cussion of the competing spin-dependent processes is the 
thickness-dependent difference between r m j n for 6 and e, 
which is qualitatively reproduced by our simulation. The 
Co-Cu interface region, which is preferentially probed by 
e, reaches the minimum magnetization earlier and also 
starts to recover earlier. Consequently, this confirms our 
above explanation of the initial dynamics by spin trans¬ 
port near the Co/Cu interface and the later dynamics 
by spin-flip scattering at the surface. We qualitatively 
reproduce the transient behavior of our MO observables 
already with such a relatively simple model up to delays 
of 400 fs when the electronic system is in non-equilibrium 
with respect to the lattice. In Fig. [5[a,b) we show two 
simulations of A 6 and Ae for d = 6, 20 nm. For 20 nm 
we obtain an initially faster (< 100 fs) reduction of e 
compared to 9 , and later at 200 fs the crossing of both 
curves. For 6 nm our simulation does not exhibit any 
difference for the two magneto-optical observables which 
is in agreement with the experimental result as far as 
for these small thicknesses the minima in Ae(t)/eo and 
A9(t)/9o are not discernible within the error bars. We 
will come back to this point in the discussion below. 

Having explained the origin of the different transient 
behavior of the two MO observables by the two separate, 
elementary spin-dependent processes of spin transport 
from the Co film into the substrate and local spin-flip 
scattering within the Co film, we simulated the spatio- 
temporal dynamics of the magnetization change Am as 
a function of position within the Co film and time de¬ 
lay. The result is shown in Fig. j5]c) and highlights the 
separate locations of the spin-flip and the spin transport 
processes in the near surface region at delays just above 
200 fs and near the Co-Cu interface earlier than 100 fs. 

Note that we used a reduced fluence in the calcula¬ 
tion compared to the experiment in order to adjust the 
calculated magnitude of demagnetization to the experi¬ 
mentally observed one. The relative difference between 
the fluence in the calculation and in the experiment is 
roughly two times which is attributed to deviations of 
the optical constants in the epitaxial film-substrate sys¬ 
tem, limitations in the determination of the laser focus 
on the sample surface in the vacuum chamber, and the 
employed simplifications in the model. 


III. DISCUSSION 

The calculated behavior with spin-flip and spin cur¬ 
rent contributions agrees with the experimental observa¬ 
tion, c.f. Figs. 12131 For t > t*, the pump induced varia¬ 
tion for 9 is stronger than for e, also in agreement with 
Fig. [2] and the dynamics are dominated by spin-flip scat- 



© 0 0.2 0 0.2 

delay [ps] 



delay [ps] 

FIG. 5. (a,b) Simulated magneto-optical response A8(t)/9o 
(solid) and Ae(t)/eo (dashed) for d = 6,20 nm. (c) Simu¬ 
lated spatio-temporal variation of the relative magnetization 
change Am for a 20 nm thick Co film on Cu(001) in a false 
color representations as a function of time delay and position 
within the film. 


tering. This behavior is only obtained if the SD contri¬ 
bution is included, the M3TM alone is not sufficient. We 
conclude furthermore that the contributions from both 
mechanisms are comparable and therefore both have to 
be considered. However, as we demonstrate here, they 
can be separated in the time domain. For d = 6 nm we 
lose the sensitivity to the competing processes in the sim¬ 
ulation, probably due to a more homogeneously demag¬ 
netized film. After all, Fig.QJa) indicates a loss of m due 
to SD at the inner interface as well as spin-flip scattering 
at the surface both on a length scale comparable to 6 nm. 
In this thin film limit the magnetization gradients result¬ 
ing from spin-flip scattering and spin transport in the Co 
film have essentially receded. They become more simi¬ 
lar to each other because the film thickness approaches 
the spin-dependent mean free path, as spins can now be 
transported into Cu also from regions closer to the Co 
surface. 

The deviation of the simulation and the experimental 
data in Fig. [3] is potentially a result of uncertainties in 
the optical constants. In addition, spin -dep endent back- 
diffusion of electrons from Cu into Co [ll|, |lj| , which is 
not included in the simulation, can enhance the demagne- 
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tization near the Co/Cu interface found in our SD simula¬ 
tions further. Spin-dependent scattering at defects could 
also contribute to the spin dynamics and would be de¬ 
termined by the defect density at the interface. Here, we 
investigate an epitaxial film on a single crystal substrate 
which is the best way to avoid such defects. As the results 
of our simulation only show a qualitative agreement with 
our experimental results, see Fig. [3l such approximations 
might be one of the reasons for this. For a quantitative 
description energy and spin dependent transmission co¬ 
efficients at the interface need to be considered, as well as 
many body renormalizations. This is, however, beyond 
the scope of this article. Likely, these aspects need to be 
taken into account for the deviation of the fluence used 
in the simulation and in the experiment. 

A punrp-induced variation of the optical constants due 
to the hot electron distribution could influence our re¬ 
sults, but we are convinced that m(z,t) dominates the 
dynamics. We argue as follows, (i) For films d > 10 nm 
the difference between 9 and e up to t = 200 fs is com¬ 
parable to later times, but with opposite sign (Fig. [2J. 
To explain this by a time dependent change of optical 
constants, a change in sign at a remaining absolute value 
would be required. This is rather unlikely since such ef¬ 
fects would decay monotonously with the hot electron 
distribution within ~ 1 ps [lg, i2al ■ (ii) With decreas¬ 
ing d the difference between 9 and e at t < t* shrinks 
faster than for t > t* (Fig. [2]), while the hot electron dis¬ 
tribution remains for all d. In contrast, spin polarized 
transport is strongly affected by the sample thickness 
[Zi EO, [3 • (iii) The difference between 9 and e remains 
up to 25 ps, see Fig. U which is too long for state fill¬ 
ing effects due to a hot electron distribution, considered 
previously [l5j]. A magnetization profile m(z,t), which is 
detected through the different ( for 9 and e, can, how¬ 
ever, remain for such a long time due to a spatial gra¬ 
dient in the lattice temperature. In fact, the gradient 
in the magnetization might persist even longer than it 
takes the lattice temperature to homogenize, as at longer 
timescales, when the excitation of the electronic system 
has already relaxed, changes in the magnetization due to 
changes in the lattice temperature are mediated by spin- 
lattice coupling, which is not included in the M3TM and 
has characteristic timescales of several ps to tens of ps. 


IV. CONCLUSION AND OUTLOOK 

We showed that the different depth sensitivity of the 
magneto-optical Kerr rotation and ellipticity can be used 
to identify spatial magnetization profiles on ultrafast 
time scales resulting from spin polarized transport or 
spin-flip scattering. We found that the laser-excited spin 
dynamics in Co/Cu(001) films are dominated by spin 
transport effects on times up to ~ 100 fs when the elec¬ 
tronic system has not yet thermalized, and by spin-flip 
scattering of thermalized electrons subsequently. 

Since the sensitivity function, which governs the depth 


sensitivity, is straight forward to calculate, the demon¬ 
strated method is readily applicable to further material 
systems. Here, analysis of metallic ferromagnets on an in¬ 
sulating substrate might provide a film substrate combi¬ 
nation complementary to the metal-metal case discussed 
here. We note that for the here presented methodology 
it will be essential to determine the complex sensitiv¬ 
ity function independently of other heterosystems by a 
full, thickness dependent analysis, because it is set by the 
film-substrate combination rather than by the film alone. 
Moreover, an improved experimental interface sensitivity, 
which can be obtained by using non-linear optical tech¬ 
niques, promises direct access to the spin transfer dy¬ 
namics across interfaces and might be rewarding in order 
to obtain deeper insight towards a full reconstruction of 
m(z , t) and a more quantitative understanding. 
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APPENDIX 

The group velocities as well as the densities of ex¬ 
cited carriers in bulk Co were calculated within the spin- 
polarized generalized gradient approximation (GGA) [44| 
to density functional theory as implemented in the Vi¬ 
enna Ab Initio simulation package (VASP) mil with the 
projector augmented waves basis sets [46j, |47T |. The cal¬ 
culations were performed with a hep unit cell of Co (lat¬ 
tice constants a ~ 2.51 A and c ~ 4.07A (3) using 
40 x 40 x 40 k-meshes for the Brillouin zone integration 
and a plane-wave cut-off of 335 eV. 

The excitation probability for an electron-hole pair 
formed by states | c,k,a) (electron) and | v,k,a) (hole) 
at crystal momentum k and with spin a is calculated 
based on Fermi’s golden rule according to 

P{k , <T, C, V ) = |\xx ‘ fi(Pc,k,a E v j zrJ F ex ), 

where \M^ va \ xx are the optical transition matrix el¬ 
ements for electric fields polarized in cc-direction, E v / C ^ 
conduction and valence band energies, and E ex = 1.55 eV 
is the laser excitation energy. We employed 0.01 eV Gaus¬ 
sian broadening for the evaluation of the excitation prob¬ 
abilities. 

The energy resolved density of excited carriers D is 
determined by 



FIG. 6. Calculated density of carriers excited in Co vs. their 
energy with respect to the Fermi level. 


^ E \M k \l x -S(E^-E^-E ex ).S(E^ k -w 

c,v,k 

where is the number of fc-points. Here, v and c 
are band indices referring to states below and above the 
Fermi level respectively. The results is depicted in Fig. [Gj 
Using the above equation we estimate the energy of pri¬ 
mary excited electrons, which yields 0.7 eV and 1.1 eV 
for the majority and minority electrons, respectively. The 
components of the group velocities perpendicular to the 
surface, which enter our above simulation, are calculated 
according to v cr (E a ) = h~ 1 d/E a dk z . 
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